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ABSTRACT 

Context. We study the feasibility of detection of neutral interstellar deuterium by the forthcoming NASA SMEX mission IBEX. 
Aims. Using numerical simulations, we study the absolute density and flux in Earth orbit of neutral interstellar deuterium and check 
its detectability using IBEX. 

Methods. Our simulations were performed using the Warsaw 3D time-dependent test-particle model of neutral interstellar gas in 
the inner heliosphere, which was specially adapted to the case of deuterium, and state-of-the-art models of the ionization field and 
radiation pressure. The modeling predicted the density, bulk velocity, and flux of interstellar D at different positions of the Earth 
during the solar cycle. We paid particular attention to the time interval in which IBEX observations will be performed. 
Results. Using our simulations, we predict a large enhancement of deuterium abundance in Earth orbit with respect to the abundance 
at the termination shock. The energy of the D atoms at IBEX will be within the energetic sensitivity band of its Lo instrument, apart 
from during a short time interval between September and November each year. Because of the specific observing geometry of IBEX, 
there will be one opportunity each year to search for I/S D, when Earth is close to ecliptic longitude 165°, i.e. in March. Assuming 
that the TS abundance of D is identical as in the Local Cloud, which is equal to 1.56 x 10~ 5 , and that the density of H at TS is 
0.11 cirT 2 s -1 , we estimate the expected relative flux to be approximately 0.015 cirT 2 s -1 , which corresponds to the local absolute 
flux of about 0.007 cirr 2 s -1 . The dependence of the expected flux on the phase of solar cycle is relatively weak. The flux scales 
proportionally to the density of deuterium at the termination shock and depends only weakly on the bulk velocity and temperature of 
the gas in this region. 
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1. Introduction 

The abundance of deuterium in the Local Interstellar Medium 
used to be a subject of debat e that appears to have been reso lved 
only a couple of years ago (Linsky et al. 2006; Linsk vl2~007l) . Its 
total value, which challenges the chemical evolution theory of 
the Galaxy, is approximately 23 ppm; in the gas phase, however, 
it is only 15.6 ppm. This knowledge is derived primarily from 
analyses of observations of UV lines of interstellar matter seen 
in the spectra of nearby stars, i.e. averaged over parsecs and in 
most cases over different clouds in the LISM. Technology avail- 
able today makes it feasible to measure the D density in the LIC 
by direct detection of neutral deuterium in situ in the inner he- 
liosphere either by observations of derivative populations (D + 
pickup ions), or by direct detection of the interstellar D atoms 
themselves. Hence, a reconnaissance of modifications of density 
and flux distribution of interstellar D inside the heliosphere, and 
of its abundance with respect to hydrogen, becomes a relevant 
and urgent task. 

The only study of the distribution of in terste l lar D in the 
heliosphere that we are aware of, is due to iFahri d!979l) . who 
performed the analysis in the appr oximation o f the perfec t hot 
model of the neutral interstellar gas dFahrl 1978b . IFahri d 1 9791) as- 
sumed that the normalized radiation pressure, expressed in units 
of solar gravity force (the so-called //-units), acting on D, is equal 
exactly to half of the normalized pressure acting on H. Fahr pre- 
dicted a gradual increase in the abundance of D with respect to 
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H with decreasing heliocentric distance and a persistence of the 
overdensity of D in the downwind region throughout the solar 
cycle. However, because of the lack of relevant data at this time, 
an important aspect of the problem was neglected. 

Interstellar D is subjected to identical ionization processes 
in the heliosphere as H. The only differences in factors shaping 
its distribution in the heliosphere appear to be its lower thermal 
velocity, because of the higher atomic mass, and a different ra- 
diation pressure, because both of the mass difference and the 
isotope effect that shifts its Lyman-a resonance wavelength by 
-0.0333 nm with respect to the center of the self-reversed solar 
Lyman-of line. It turns out that the latter effect has far-reaching 
consequences for the distribution of abundances, densities, and 
fluxes of interstellar D in the heliosphere, especially in the con- 
text of its potential detection by a forthcoming NASA SMEX 
mission IBEX, scheduled for laun ch in the Summer of 2008 
dMcComas et alJl2004l2005Ll2006l) . 

We begin with a brief discussion of the factors affecting the 
kinematics of neutral D atoms in the inner heliosphere and we 
present the simulations performed. Subsequently, we present the 
distribution of density, bulk velocity, flux, and abundance of neu- 
tral interstellar D along the Earth orbit and their variations during 
the solar cycle. We conclude by presenting the conditions under 
which interstellar D will be detectable by the IBEX mission. 

2. Model and simulations 

Neutral interstellar deuterium atoms in the heliosphere are colli- 
sionless and obey the equation of motion that includes the attrac- 
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Fig. 1. Compensation factors of solar gravity due to Lyman-a ra- 
diation pressure for H and D atoms for solar minimum and max- 
imum conditions as a function of radial velocity of the atoms. 
The profiles, scaled in km/s of equivalent Doppler shift and in fi 
units, are shown according to the models defined in Eq.© and 

(so- 

(solar max.). Approximate 
spectral regions relevant for D and H are marked in gray. 



Eq.©, for the net flux / tot equal to 3.53 x 10 11 cm 2 s 1 
lar min.) and 5.49 x 10 11 cm -2 s _I 



tion by solar gravity and repulsion by the solar EUV radiation in 
the Lyman-o- line: 



d 2 r/dt 2 = -GM[l -fi(v r (t) , / t ot (»))] r/\r\ : 



(1) 



where r (?) is the position vector of the atom with respect to the 
Sun at time t and G M the gravity constant times solar mass. The 
compensation factor of solar gravity due to radiation pressure 
p depends on the radial velocity of the atom v r and the line- 
int egrated solar Lyma n- a flux / to t (t), which was s hown for H 
bylTarnopolski d2007l) andlTarnopolski & Bzowski (2007) using 
obser vations of the solar Lyman-o' line profiles by lLemaire et al.l 
(120021) : 

H (v r , / to t) = A(l+B / tot ) exp (-Cv 2 ) (2) 
x [l + D exp (Fv, - Gvj) + H exp (-Pv r - Qvj)] 

For deuterium, the compensation factor fi D is derived from the 
formula for hydrogen (f2]i by changing the variable v r to v r - 
82.1201 k m/s and scaling th e result down by the D/H atomic 
mass ratio (Tarnopolski 2007) to obtain: 

\xd (v r , /tot) = a (1 + b / to t) [exp {-c v r -d v 2 ) + 

+ /exp(-gv r -/iv 2 )+pexp(-^v f -rv 2 )], (3) 

with the following parameters: 

a = 4.9469 x 10 7 , b = 4.5694 x 10~ 4 , 

c = 2.3603 x 10~ 2 , d = 3.8967 x 10~ 4 , 

/ = 5.6579 x lO- 4 , g = 0.10795, h = 3.7205 x 10~ 4 , 

p = 0.52459, q = 6.2923 x 10~ 3 , r = 3.8312 x 10~ 5 . 

The normalized radiation pressure (measured in ju-units) act- 
ing on D would be just half of the pressure acting on H if the 
solar Lyman-a line were flat. An isotope effect creates a small 
difference in the Lyman-a resonance wavelength of the H and D 
atoms. This difference of -0.0333 nm, albeit small, is responsi- 
ble for the phenomena in the heliospheric gas that we discuss in 
the paper. Since the solar Lyman-alpha line is not flat, as demon- 
strated in Fig. [1] the spectral fluxes at the resonance Lyman- 
a wavelengths of H and D atoms are different; when H and D 
atoms have identical radial velocities, the normalized radiation 



pressure force that they experience usually differs therfore by a 
factor other than 0.5. When the radial velocity of a D atom mov- 
ing towards the Sun is larger than the radial velocity of a H atom, 
the differences in radiation pressure increase still more. 

The range of radial velocities of helios pheric H atoms from 
the primary and secondary populations (Bzo wski et al.l [19971 
Izmodenov 2001) position them within the central self-reversal 
of the profile: the usual assumption that the radiation pressure 
does not depend on ra dial velocity therefore appears reasonable 
although not perfect (Tarnopolski & Bzowskill2007i) . By con- 
trast, the radial velocities of D fall on the blue slope of the so- 
lar line, so the radiation pressure is highly asymmetric with re- 
spect to the zero radial velocity. Before perihelion, the atoms are 
subject to a relatively high radiation pressure - a little higher 
than half of the radiation pressure acting on the corresponding 
H atoms. This ensures that the D atoms accelerate to a higher 
speed on the pre-perihelion leg of their orbits than the H atoms 
do. After perihelion, this speed is gradually reduced because the 
D atoms wander to the wing region of the solar line where the 
resulting compensation of solar gravity is appreciably reduced 
with respect to the values at the inbound wing, and the net force 
is therefore highly attractive. Hence the D atoms are strongly de- 
celerated at the outbound leg of their trajectories and leave the 
solar vicinity with a lower speed than they had when approach- 
ing the Sun, especially during solar minimum. For example, an 
atom that passes the Sun at a distance of 1 AU traveling at a ve- 
locity of approximately 32 km/s and had a velocity of 22 km/s 
at 100 AU in front of the Sun, has a velocity of about 15 km/s at 
100 AU behind the Sun at solar minimum. At solar maximum, 
it is caught within a bound orbit, where it becomes ionized. For 
the atoms not trapped in bound orbits, the entrance speed would 
of course equal the exit speed, if the radiation pressure force had 
been symmetric with respect to zero radial velocity. 

The simulations of densities, bulk velocities, and fluxes 
of D at 1 AU were performed using the Warsaw test-particle 
3D time-de pendent model of interstel l ar gas in the inner 



heliosphere dR ucinski & Bzowskil 119951: IBzowski et ail 119971: 
Bzo wski et alj 120081) . whic h had been specially adapted for deu- 
terium (Tarnopolski 2007). The classical approach of the hot 
model dThomaslll978l) was used, but instead of Keplerian mo- 
tion the atoms were assumed to follow trajectories modeled us- 
ing numerical solutions of the equation of motion Eq.(Q]) in ad- 
dition to Eq.Q. The v ariations in the charge-exchange rate with 
both latitude and time (IBzowski et alj200a : IBzow ski 2001), and 
the electron impact ionization in a stat ic, spherically symmet- 
ric approximation (Bzowski et al. 2008), were all taken into ac- 
count. The radiation-pressure factor fi was governed by the net 
flux / tot in the solar Lyman-q line, whic h was taken from the 
SOLAR 2000 model (Tobiska et a j]|2000|) and approximated by 
the analytic formula f rom lBzowskil (1200 ll) with parameters from 
IBzowski etai] (2008). Details will be shown in a future paper. 

The simulations were carried out for the ~ 1 1-year cycle of 
solar activity 1986-1998 and followed the Earth position during 
its yearly motion for the 5-th day of each month; the calcula- 
tion point in June and December of each year therefore coin- 
cided with the projection of the inflow direction on the Earth 
orbit. The upw ind d i rectio n was assumed to equal the direction 
calculated by Witte (2004), using in situ observations of inter- 
stellar He. The density of interstellar H at the termination s hock 
(TS) was assumed to equal 0.11 cm -3 (Bzowski et al. 20081) and 
the abundance of D with respect to H at TS was assumed to 
be identical to that of th e gas phase of the Local Interstellar 
Cloudy = 1.56 x 10~ 5 (iLinskvet al.ll2006l) . The bulk veloc- 
ity and temperature of the gas at TS were equal to 22 km/s 
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D density & D/H ratio at Earth orbit 
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Fig. 2. Density of interstellar D along the Earth orbit during the 
solar cycle (black, left-hand scale) and abundance of D with re- 
spect to H (gray, right-hand scale) for the D/H TS abundance 
1.56 x 1(T 5 and H density 0.11 cm 3 . 



and 12000 K, respectively, close to the average values of these 
parameters from the primary and secondary populations of in- 
terstellar hydrogen from the Moscow MC model adopted by 
iBzowski et al.l d2008l) to obtain the H density used in the sim- 
ulations, andjnagreement with values inferred for interstellar H 
at TS bv lCosta et all (Il999h . 

The problem of filtration o f D thr ough the interface has not 
yet been addressed. Fah r et all ( 1 19951) simulated the filtration of 
interstellar O assuming that it is almost in charge exchange res- 
onance with protons, as is the case for D; they discovered that 
the balance between the primary and secondary populations at 
TS is different from the case of hydrogen, which produces net O 
parameters at TS that are closer to th e parameters in the unper- 
turbed gas. lOsterbart & Fahrl ( 1 1992b showed a breakdown be- 
tween H populations that originate in various collision hierar- 
chies which, together with the oxygen study, suggest that in the 
case of D the secondary population may have different parame- 
ters from the secondary population of H. These conclusions were 
supp orted by simula t ions o f O transmission through the interface 
by llzmodenov et all dl997l) . who found a transmission factor of 
O equal to 0.7. To assess the robustness of our results, we there- 
fore completed a parameter study of the 1 AU density of D as a 
function of the bulk velocity and temperature at TS. 

With deuterium simulations, corresponding simulations of 
hydrogen were also performed. 

3. Results 

Our simulations indicate that deuterium abundance in Earth or- 
bit is elevated at least twofold with respect to the abundance at 
TS and varies by a few orders of magnitude along the Earth or- 
bit from upwind to downwind, where it reaches its highest value 
(see Fig.[2j). It also varies during the solar cycle. The explanation 
of these reults is a strong v ariability of neutral in terstellar hydro- 
gen density in Earth orbit (IBzowski et al.| [l996). The density of 
deuterium has a relatively mild variability (at upwind, by a fac- 
tor of 0.5 from solar minimum to maximum), as shown in Fig. [2] 
The bulk speed of deuterium varies weakly during the solar cycle 
and remains almost constant along the Earth orbit (30 - 32 km/s) 
apart from a tail region (gray line in Fig. [3}. Consequently, the 
net flux of interstellar D shows an annual amplitude that is only 
a little larger than the amplitude of density (see the gray line in 
FiglU). The amplitude decreases with the increase of solar activ- 
ity, but the variabilities in the flux in the upwind and crosswind 
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Fig. 3. Bulk speed of interstellar D at Earth orbit during the phase 
of solar cycle corresponding to the planned operations of IBEX 
(gray) and relative speed of D with respect to Earth (black). 
Vertical lines mark the interval corresponding to IBEX obser- 
vations, whereas the gray bands mark the intervals when the D 
beam should be visible to IBEX. 



D flux at Earth orbit, D/H= 1 .56x 1 0~ 5 , n H . T s =0.11 cm" 3 
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Fig. 4. Net flux of neutral interstellar D expected at Earth orbit 
(black) and the corresponding flux relative to Earth-bound ob- 
server (gray) as function of time. Vertical lines mark the interval 
corresponding to the planned IBEX observations, and the gray 
bands mark the intervals when the D beam should be visible by 
IBEX. 



portions of Earth orbit are relatively small: the solar max/solar 
min flux ratios are equal to ~ 0.4 and ~ 0.25, respectively. 

The time profile of the velocity of atoms hitting a detector on 
an Earth-bound satellite is appreciably modified by the proper 
velocity of the spacecraft, which is assumed here to equal the 
Earth velocity, as shown by the black line in Fig. [3] Since the 
detection efficiency usually increases with energy, the optimal 
time for detection should be an interval close to the 40th day 
of each year, when the relative velocity is the highest. From the 
viewpoint of flux strength, the optimal time is about 135-th day 
of each year (black line in Fig. |4j. 

IBEX will be a spin-stabilized Earth satellite that is designed 
to observe energetic neutral atoms (ENA) in the 0.01 - 5.9 keV 
energy band, with the centers of relevant energy bands at 0.015, 
0.029, 0.056, and 0.107 keV and sensors looking perpendicu- 
larly to the spin axis directed at the Sun. It is planned to operate 
for at least two years. T he observation techn ique for the slowest 
atoms was discussed by Mob ius et al.l (12001). The opening angle 
of the IBEX-Lo instrument, best suited for observations of inter- 
stellar D, will be 7° and the rotation axis of the spacecraft will 
be repositioned approximately every seven days to maintain its 
offset from the Sun at no more than ~ 4°. Thus the beam of in- 



4 



Tarnopolski & Bzowski: Detectability of interstellar D by IBEX 



D relative energy ± thermal spread at Earth orbit 
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Fig. 5. Relative energy of deuterium gas with respect to the mov- 
ing Earth + thermal spread. Horizontal lines mark the centers of 
the lowest IBEX energy channels (except the lowermost one, 
which marks the lower sensitivity limit), the vertical lines mark 
the interval corresponding to the planned IBEX observations. 
The "sweet points" for deuterium detection are marked by the 
gray vertical bands. 

terstellar D atoms will be observable by IBEX when the ecliptic 
longitude of its vector of relative velocity with respect to Earth is 
inclined by 90° ±7° with respect to the Sun-Earth line. Hence, an 
opportunity to detect interstellar D atoms can occur potentially 
only twice during each year: once when IBEX is between ~ 105° 
and ~ 165° ecliptic longitude (between January and March), and 
the other when the spacecraft is located symmetrically with re- 
spect to the upwind longitude 254.68°, i.e. between October and 
December. But during the latter season the relative speed be- 
tween the D flow and IBEX will be so slow that the atoms will 
be unable to exceed the lower boundary of the energetic sensitiv- 
ity band; only the January-March opportunity is therefore viable, 
when the Earth travels against the gas flow (Pig. |5j. 

The beam of D atoms has a thermal spread in energy whose 
magnitude varies with position along the Earth orbit. The ther- 
mal spread projected onto an arbitrary line, defined by the unit 
vector e(r), is calculated from the formula: 

where u is the local bulk velocity of the beam. Fig. [5] shows the 
relative energy of the D atoms incoming to IBEX within uj (r) 
from the local velocity u (r). These atoms are located within the 
IBEX sensitivity limit throughout the year apart from a short 
interval between August and November. The beam portion con- 
tained within one thermal velocity uj from the local bulk speed 
u will be within the IBEX sensitivity range precisely during the 
first half of the calendar year. The entire energy range of the 
beam will drop below the detectability limit of IBEX-Lo about 
the 230th day of each year and reemerge about DOY 320. 

The lower limit to the flux detectable by the IBEX-Lo instru- 
ment is conservatively estimated to be 0.1 - 1.0 at. cm _2 s _1 sr _1 
[Mobius, 2008, private comm.]. During the January-March ob- 
servations opportunity, the net flux observed by IBEX is ex- 
pected to increase from 0.002 cirT 2 s _I to 0.015 ctrT 2 s _1 , 
which corresponds to the absolute fluxes equal to 0.001- 
0.007 cirT 2 s _1 . Under the assumptions adopted in this study, 
the expected flux of D is therefore at the detection threshold. 

Calculations suggest that this result depends weakly on pos- 
sible modifications of temperature and bulk velocity of neu- 
tral D at the termination shock. For solar minimum conditions, 



a change in the gas temperature at the termination shock by 
+6000 K (i.e., by +27%) from the assumed value of 12 000 K 
produces a change in the density at 1 AU crosswind by +25%, 
and a change in bulk velocity at the termination shock of +4 km/s 
(i.e., by +18%) from 22 km/s results in a change in the 1 AU 
density by ±10%. For solar maximum conditions, these varia- 
tions will be approximately twice as large when the temperature 
is varied, and will remain unchanged when the bulk velocity is 
varied. An increase in the temperature and/orbulk velocity at TS 
leads to an increase in deuterium density at 1 AU. On the other 
hand, since the local density and local flux are directly propor- 
tional to the density at the termination shock, then if the filtration 
of D within the heliospheric interface is weaker than in the case 
of hydrogen, or if the density of D in the LIC is higher than in- 
ferred from the astrophysical observations, then the density at TS 
might be higher and the flux at 1 AU could exceed more readily 
the IBEX-Lo detection threshold. 

4. Conclusions 

We performed extensive simulations of the density, velocity, and 
flux of neutral interstellar deuterium at Earth orbit in an attempt 
to check its detectability by IBEX. The conclusions from this 
study are the following. 

1 . The absolute density of neutral interstellar D at Earth orbit 
changes during the solar cycle more weakly the correspond- 
ing H density and its abundance increases appreciably with 
respect to the abundance at the termination shock. 

2. The bulk velocity of D at Earth orbit depends weakly on the 
ecliptic longitude (apart from a narrow cone in the down- 
wind region), so variations of the absolute flux follow the 
variations of density. 

3. The relative energy of D, and its thermal spread, are within 
the IBEX sensitivity band apart from an interval between 
August and November each year, but because of the mea- 
surement geometry the only observation season during the 
year occurs when the Earth is traveling "against the wind" of 
interstellar D at Earth orbit, i.e. between January and March. 

4. The expected relative flux is equal to 0.015 cirT 2 s , which 
corresponds to the absolute flux equal to 0.007 cirT 2 s _1 , as- 
suming the H density at TS equal to 0.11 cirT 3 and a D/H 
abundance in the gas phase at TS equal to 15.6 ppm. 

5. Under these assumptions, the expected flux of D appears to 
be at the sensitivity limit of the IBEX-Lo instrument. 
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